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Abstract 

We report on the search for the factorization suppressed decays B — > XcO ^ and B — > Xc2 , 
with XcO and Xc2 decaying into J/tp 7. We use a sample of 124 million BB events collected with 
the BABAR detector at the PEP-II storage ring at the Stanford Linear Accelerator Center. No 
significant signal is found and upper bounds for the branching fractions are obtained. All results 
are preliminary. 
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1 INTRODUCTION 



Hadronic decays of heavy mesons are not precisely described, despite the electroweak nature of the 
quark decay, because the initial and final states consist of mesons, not of quarks. The factorization 
scheme allows one to make some predictions though. Factorization assumes that a weak decay 
matrix element can be described as the product of two independent hadronic currents. Under 
the factorization hypothesis, B —> c K^*> decays are allowed when c = J/ip , ip(2S) or Xci, but 
suppressed when c = XcO or Xc2 CO- I* 1 lowest order heavy quark effective theory, there is no J > 2 
operator to create the tensor Xc2 from the vacuum. The decay rate to XcO is zero due to charge 
conjugation invariance 

Belle has recently observed B — » XcO K + decays, with XcO - » k + ir~ or K + K~, with a 
branching fraction surprisingly large based on the expectation from factorization and measurements 
of the Xci branching fraction. BaBar has confirmed the observation A a with a branching fraction 
somewhat lower than, but compatible with, that measured by Belle. 

In this document we attempt the detection of B — > Xc,i K^*\ i = 0, 2, using the radiative Xc —> 
J/ip 7 decays. 

2 THE BABAR DETECTOR AND DATASET 

The data used in this analysis were collected with the BABAR detector at the PEP-II storage ring. 
They represent an integrated luminosity of 112.4 fb _1 of data taken at the T(4S) resonance. 

The BABAR detector is described elsewhere (Hj- Charged particles are detected with a five-layer, 
double-sided silicon vertex tracker (SVT) and a 40-layer drift chamber (DCH) with a helium-based 
gas mixture, placed in a 1.5-T solenoidal field produced by a superconducting magnet. The charged- 
particle momentum resolution is approximately {5pt/pt) 2 = (0.0013pr) 2 + (0.0045) 2 , where px 
is the transverse momentum in GeV/c. The SVT, with a typical single-hit resolution of 10 /im, 
measures the impact parameters of charged-particle tracks in both the plane transverse to the 
beam direction and along the beam. Charged-particle types are identified from the ionization 
energy loss (dE/dx) measured in the DCH and SVT, and from the Cherenkov radiation detected 
in a ring-imaging Cherenkov device (DIRC). Photons are identified by a CsI(Tl) electromagnetic 
calorimeter (EMC) with an energy resolution a{E)/E = 0.023 • [Ej GeV)" 1 / 4 0.019. Muons are 
identified in the instrumented flux return (IFR), composed of resistive plate chambers and layers 
of iron, which return the magnetic flux of the solenoid. 

3 ANALYSIS METHOD 

The channels considered here are B — > Xc with Xc — > J/i/> 7> Jfy ~~ * e + e _ or K is K + 

or K° (vr+ tT); K*° -» K+ tT or K®, vr°; K* + -> K + vr° or K° vr+; and tt° -»• 7 7. 

Multihadron events are selected by demanding a minimum of three reconstructed charged tracks 
in the polar-angle range 0.41 < 0^ < 2.54 rad. Charged tracks have to be reconstructed in 
the DCH and are required to originate at the beamspot: within 1.5 cm in the plane transverse 
to the beam and 10 cm along the beam. Events are required to have a primary vertex within 
0.5 cm of the average position of the interaction point in the plane transverse to the beamline and 
within 6 cm longitudinally. Electromagnetic depositions in the calorimeter in the polar-angle range 
0.410 < Oiab < 2.409 rad that are not associated with charged tracks, have an energy larger than 
30MeV, and a shower shape consistent with a photon are taken as photon candidates. A total 
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energy larger than 4.5 GeV in the fiducial regions for charged tracks plus neutrals is required. To 
reduce continuum background, we require the normalized second Fox- Wolfram moment R2 (HI of 
the event, calculated with both charged tracks and neutral clusters, to be less than 0.5. Charged 
tracks are required to be in polar-angle regions for which the PID efficiency is well-measured. For 
electrons, muons, and kaons the acceptable ranges are 0.40 to 2.40, 0.30 to 2.70 and 0.45 to 2.50 
rad, respectively. 

Event selection was optimized by maximizing e / y/B, where e is the signal efficiency, and B the 
number of events, after selection was applied. Candidates for J/tjj mesons are reconstructed in the 
e + e~ and decay modes, from a pair of identified leptons that form a good vertex. Muon 

candidates are identified using a neural network selector and are required to pass a loose selection 
for the first muon candidate and a very loose selection for the second muon candidate. Electron 
candidates are selected using a likelihood selector and are required to pass a loose selection. For J/tp 
— > e + e~ decays, electron candidates are combined with photon candidates in order to recover some 
of the energy lost through bremsstrahlung. Photons are required to be within 35mrad in polar 
angle from the electron track, and to have an azimuthal angle intermediate between the initial track 
direction (estimated by subtracting 50mrad opposite to the bend direction of the reconstructed 
track) and the centroid of the EMC cluster arising from the track. 

The lepton-pair invariant mass has to be between 2.95 and 3.18 GeV/ c 2 for both lepton flavors. 
The remaining background is mainly due to genuine J/ip 's. Candidates for K® consist of oppositely- 
charged tracks with invariant mass between 487 and 510 MeV/c 2 and are required to satisfy vertexing 
conditions. The K® flight length has to be greater than 1 mm, and its direction must form an angle 
with the K® momentum vector in the plane perpendicular to the beam line that is less than 0.2 rad. 
Charged kaon candidates are identified using a likelihood selector and are required to pass a tight 
selection. Photon candidates as defined above are used also for the reconstruction of 7r° — ► 77 
decays. A 7r° candidate consists of a pair of photon candidates with invariant mass in the interval 
117 - 152MeV/c 2 and momentum larger than 350 MeV/c. K* candidates must have a Ktt invariant 
mass in the range 0.85 - 0.94 GeV/ c 2 around the nominal K* (892) mass [7j. The J/ift, K®, and 7r° 
candidates are constrained to their corresponding nominal masses [J]. 

The Xc candidates are formed from J/ip and photon candidates. The photon is required to fulfill 
the same shower shape requirement mentioned above, have an energy larger than 0.15 GeV, and 
not be part of 7r° candidates in the mass range 0.125 - 0.140 GeV/c 2 . 

The Xc and candidates are combined to form B candidates. Two kinematic variables 

are used to further remove incorrect B candidates. The first is the difference AE = E* B — E£ eam 
between the -B-candidate energy and the beam energy in the T(4<S) rest frame. In the absence 
of experimental effects, reconstructed signal candidates have AE = 0. The second is the beam- 
energy-substituted mass 771ES = (Eteam ~P*b) 1 ^ 2 - The energy substituted mass tties should peak 
at the B meson mass 5.279 GeV/c 2 . 

For the signal region, AE is required to be between — 35MeV and +20MeV for channels involv- 
ing a 7T°, and to be between —20 MeV and +20 MeV otherwise. If several B candidates are found in 
an event, the one having the smallest \AE\ is retained, toes is required to be in the 5.274 - 5.284 
GeV/c 2 range. 

Studies using simulated samples show that most of the background events in the Xc K* channels 
are due to non-resonant (NR) B — ► Xc (J/ifi l)K ir decays. Also the observation of the suppressed 
Xc2 could be complicated by the presence of the prominent Xci peak. Therefore the search is 
performed by the observation of the spectrum of the mass difference m^+^- 7 — mp+g- . It was found 
from Monte Carlo simulation that after the non-resonant events were removed from the sample, 
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the expected number of genuine Xc ~^ J/ip 7 decays was extremely small, 0.2 ± 0.2 for the Xc2 K*° 
(K + ir~) and Xc2 K* + (K + ir°) modes, and 0.0 ± 0.2 for all the other Xc2 modes and all the XcO 
modes. 

The efficiencies obtained from fits to the mass difference distribution for exclusive samples, 
where one B decays to the final state under consideration and the other inclusively, are given in 
Table H Note that Xd exclusive simulated samples are used in the place of Xc2 , that were not 
available. 

The Xc2 has a negligible natural width and is therefore fitted with a Gaussian. The XcO has 
a natural width T = 10.1 MeV comparable with the measurement resolution a ~ 10MeV/c 2 , and 
therefore the XcO peak is fitted with a Voigtian, the convolution of a Breit-Wigner and a Gaussian 
distribution. 



Table 1: Efficiencies from fits of the distribution of m^+£- 7 — rri£+£- for exclusive samples. 





Xcl 


XcO 


K* u (K+ O 

K*° (K° s 7T°) 
K*+ (K+ 7T°) 
K*+ (K° s 7T+) 

K+ 


0.071 ± 0.001 
0.031 ± 0.001 
0.036 ± 0.001 
0.065 ± 0.001 
0.144 ± 0.001 
0.158 ± 0.001 


0.066 ± 0.001 
0.010 ± 0.000 
0.031 ± 0.001 
0.062 ± 0.001 
0.117 ± 0.002 
0.126 ± 0.001 



We corrected for the presence of non-resonant decays under the K* peak in the following way: 
the mi+i-^ — mg+i- distribution for the events on the plateau; i.e. 1.1 < niK-w < 1.3GeV/c 2 , is 
subtracted from the m£+^- 7 — mg+g- distribution for the events in the signal region 0.85 < < 
0.94GeV/c 2 , after rescaling by a factor r = 0.26 ± 0.04, where r is the ratio obtained from Monte 
Carlo simulation of non-resonant events under the peak compared to the plateau. The branching 
fractions were then computed from: 



N 

BF= (1) 

N B x e x / K ' 

where N is the number of events obtained from fitting the m^+ f - 7 — m e + e ~ distribution, is 
the number of BB events, e is the selection efficiency and / is the secondary branching fraction of 
the B daughters. Examples of fits to the "generic" BB Monte Carlo (MC) sample, that contains 
a simulation of inclusive T(4S) — > BB decays, can be seen in Fig. ^ 

The free parameters in the fits are: a linear background, the overall mass difference scale, the 
resolutions of the gaussian taken to be the same for the 3 x' s an d the amplitudes of the peaks. The 
fixed parameters are the natural width of the XcO an d 2 mass differences, all taken from PDG. 

With such fits, it was checked that the non-resonant events were subtracted correctly, and that 
the proximity of the Xcl was not inducing any visible bias on the measurement of the nearby Xc2- 



4 SYSTEMATIC STUDIES 

This measurement is affected by the following set of systematic uncertainties: 
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Figure 1: Distribution of m i + i - 1 — mi+g- (GeV/c 2 ) for generic MC samples. 



Overall uncertainty on the number of B events, 1.1%. 

Uncertainy on the secondary branching fractions: from PDG [7] (dominated by the relative 
uncertainty of the branching fraction of the radiative decay of the x, 11.9 and 8.5% for the 
Xco and Xc2, respectively). 

Tracking: 1.3% per track. 

K^: a 2.5% uncertainty. 

Neutrals: 2.5% per "x" photon, 5.0% per tt°. 

An overall 3% uncertainty on particle identification correction. 

Selection cuts: For each mass peak and for AE, the uncertainty of the MC-to-data shift in 
central value and in width are measured on the well populated Xci channels and are used 
to vary the selection cuts, by 1 a. The corresponding efficiency variation, estimated on the 
exclusive sample, is the induced contribution to the systematics. The central value and width 
induced systematics are estimated independently, and are added quadratically below. 

The results for Xci MC sample and Xco MC sample are quite close to each other; an average 
value is used for both. 

The ratio of B° to B + production in T(4S) decays is assumed to be unity. The related 
uncertainty is small [H] and is neglected here. 
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• The NR component is probably in an S-wave Kir state, as was observed in the J/ip K* system 
jSj, with an unknown relative phase <f> wrt the main K* (892) P-wave peak. 

It is possible that no signal is found in the channels under consideration in this section. 
Therefore the systematics due to the unknown relative phase is here estimated with a MC- 
based method. 

The K-7T invariant mass is fitted with an amplitude that is the sum of a non-relativistic Breit- 
Wigner and a real amplitude that corresponds to a polynomial (parabolic) distribution for 
the NR (Fig. EJ). 

2 

(2) 

where a and b are real quantities. The slow variation of the phase of the S wave with rriKn is 
neglected here. 



p(m K 7r) 



mpdg - m K7T - iT/2 



600 



400 



200 







\ 








1 

1 
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V 




L 



mK7t (GeV/c ) 

Figure 2: The distribution of rriKn for the generic MC sample. A fit with a coherent sum of a 
polynomial NR and a non-relativistic Breit-Wigner is overlaid. 



The free parameters in the fit are the 3 degrees of freedom of the NR parabola, the magnitude 
of the signal, and the relative phase (p. As the high mass plateau is dominated by the NR 
contribution, no attempt is made to subtract the few combinatorial events. The fact that eft 
is unknown is dealt with by randomly generating samples of events distributed as above for 
each value of <j>, and applying the NR subtraction as described above. The number of events 
N thus measured is normalized to that generated with the value (f>o of <f> obtained in the fit, 
and the ratio R = N \<p) / 'N (4>o) is plotted as a function of 4> in Fig. |3J The medium value is 
1.44 with maximal relative extention ± 35%, giving an RMS relative uncertainty of ± 20%. 

• In the case of Xc2, the efficiency depends on the intensity fractions to various polarization 
states, due to the variation of the detection efficiency with the angles describing the decay. 
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Figure 3: Number of events measured, after NR subtraction, in Toy MC samples, as a function of 
<fi, normalized to the number obtained with the phase fitted on the generic sample (shown by the 
vertical line). 
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Figure 4: Signal efficiency as a function of the helicity angles for K + ir channel. \ is the angle 
between the decay planes of the \c and of the K* . 
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The efficiency is mainly sensitive to the K* decay helicity angle, (Fig. |1J) due to soft pions 
for small values of 9 k* ■ 

The selection efficiency therefore depends, to 1st order, on the polarization of the K* popu- 
lation, throught the angular distribution: 



1 dr 



T d cos 9k* 



(1 - cos 2 9 K *) + l^ol^cos 2 ^* - 1) 



(3) 



where |^4o 1 2 describes the fraction of longitudinal K* polarization. The efficiency is: 

1 dr 



r dcos 9k 



-e(9K*)d cos 9k* = a + \Aq\ b, 



(4) 



where 



3 f 

a = - I (1 — cos 2 9k*)£(9k*) sin 9k*o\9k* 



(5) 



b=-f (3cos 2 # K * - l)e{9 K *)sm9K*d9K*- (6) 

The values of a and b are obtained from the two above equations and from the parametrisation 
s(9k*) extracted from Figure |I] and are shown in Table |2J In the case no signal is observed, 
the polarization is unknown, and we estimate the efficiency as (a + 0.56) ± (|6|/\/l2)- 



Table 2: Coefficients for the calculation of amplitude dependent average efficiency for the Xc2 K* 
channels (%). 





a b Efficiency Fract. uncert. 


K*" (K+ 7T") 

K*° (K° s 7T°) 

K*+ (K+ 7T°) 

K* + (K° s 7T+) 


8.68 -1.40 7.98± 0.40 5.1 
4.25 -1.66 3.43± 0.48 14.0 
5.05 -1.79 4.16± 0.52 12.4 
7.83 -1.84 6.92± 0.53 7.7 



As usual the effect is stronger on channels having a tt in the final state, as the larger 
background results in harsher cuts during the optimization process. 

A summary of the multiplicative contributions to the systematics can be found in Table El 
In addition to these multiplicative contributions comes a contribution from NR background 
subtraction. The contribution of the uncertainty of r is given in Table 01 

5 PHYSICS RESULTS 

Fits on the data clearly show the presence of the factorization allowed Xci, but no signal, within 
uncertainty, for the factorization suppressed XcO and Xc2 (Table El Fig. EJ. 

Non-resonant events subtraction has been applied. The phase-related systematics estimated on 
the MC sample are used for the data. 
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Table 3: Systematics: summary of the multiplicative contributions: relative uncertainties (%). 





K*° (K+ vr") 




(K+ TT°) 


V- I1 s " ) 


K + 




T\ poiiirhinp" 

J — ' \j\J LLJ.A Li J.J.J. tL. 


1.1 


1.1 


1.1 


1.1 


1.1 


1.1 


Tracking 


5.2 


2.6 


3.9 


3.9 


3.9 


2.6 


6 




2.5 




2.5 




2.5 


Neutrals 


2.5 


7.5 


7.5 


2.5 


2.5 


2.5 


Particle identification 


3.0 


3.0 


3.0 


3.0 


3.0 


3.0 


Cut variation () 


1.0 


1.5 


1.3 


0.8 


0.6 


0.5 


Cut variation width 


7.6 


13. 


11.5 


8.2 


6.5 


6.3 


MC stat 


1.4 


2.9 


1.7 


1.8 


1.3 


1.3 


phase 


20.0 


20.0 


20.0 


20.0 






XcO Sec. BF 


11.9 


11.9 


11.9 


11.9 


11.9 


11.9 


Total for XcO 


25.4 


28.3 


27.6 


25.5 


14.8 


14.6 


Xc2 Sec. BF 
Polar 


8.5 
5.1 


8.5 
14.0 


8.5 
12.4 


8.5 
7.7 


8.5 


8.5 


Total for Xc2 


24.5 


30.5 


29.1 


25.3 


12.2 


12.0 



Table 4: Systematics on the measured BF's on the generic samples due to NR subtraction, (in 
units of 10- 4 ). 





Xc2 XcO 


(K+ 7T-) 

(K° 7T°) 
(K+ 7T°) 

(K° s vr + ) 


0.0 0.5 
0.2 1.3 
0.1 1.4 
0.1 1.7 



Table 5: Number of events from fits of the distribution of m^+^ 7 — rri£+£- for the data. (NR 
subtracted). 





Xc2 XcO 


X* U (if + 7T-) 
if* (*T° 7T°) 
(K+ 7T°) 
(JfO 7T+) 
K+ 

K° 


0.3 ± 1.1 2.1 ± 2.5 
-1.7 ±2.0 1.0 ±0.9 
-1.8 ± 0.6 0.1 ± 2.9 
-0.2 ± 1.2 12.3 ± 3.7 
6.4 ± 4.8 15.1 ± 7.6 
2.8 ± 2.6 4.5 ± 4.0 
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Figure 5: Distribution of m f + f -^ — m f + f - for data. Top: raw data; Bottom: NR subtracted. 
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Combining the measurements of the K* sub- modes, and under the reasonable approximation 
that the multiplicative efficiencies between each K* sub-mode are fully correlated, we obtain the 
branching fractions for the suppressed modes listed in Table |H1 

The results for the allowed Xd are found to be compatible with those of an analysis 
optimized to the relevant BF, in contrast with this one. 

Table 6: Measured Branching fractions (in units of 10 ). 





Xc2 


XcO 


K *o 


0.02 ± 0.07 ± 0.10 


3.09 ± 3.35 ± 1.27 


K *+ 


-0.37 ± 0.15 ± 0.20 


27.0 ± 11.2 ± 9.0 


K+ 


0.15 ± 0.11 ± 0.12 


7.49 ± 3.76 ± 1.09 


K° 


0.18 ± 0.16 ± 0.12 


5.96 ± 5.39 ± 0.88 



Upper bounds on the BF's, at 90% confidence level (CL) are obtained using a simulation, 
assuming gaussian statistics for the statistical uncertainties and taking into account the systematic 
uncertainties (Table Q). It has been assumed that the BF can only be positive. 

Table 7: Branching fractions: upper bounds at 90% CL. (in units of 10~ 4 ). 





Xc2 


XcO 




0.22 


8. 


K *+ 


0.14 


45. 


K+ 


0.36 


12. 


K° 


0.44 


13. 



6 SUMMARY 

The upper limits obtained for decays to Xc2 are more than one order of magnitude lower than the 
branching fractions of the factorization allowed decays and of the already observed B — * 

XcO K+ 

decays. All results are preliminary. 
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